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ABSTRACT 


Ultraviolet  emissioixs  from  the  earth’s  ionosphere  in  the  wavelength 
range  1850  A  to  2070  A  are  analyzed.  Specifically,  thirteen  Lyman-Birge- 
Hopfield  (LBH)  band  emissions  from  molecular  nitrogen  (N2)  and  the  S/ 
and  e  bands  of  nitric  oxide  (NO)  are  modelled  with  synthetic  spectra  to 
determine  the  atmospheric  density  of  NO.  Corrections  to  the  Franck-Condon 
factors  for  the  LBH  emissions  and  the  NO  emissions  are  derived.  Column 
emission  rates  are  determined  for  N2  and  compared  with  computer  models. 
Vacuum  calibration  of  the  MUSTANG  ultraviolet  spectrograph  is  presented. 
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I.  INTRODUCTION 


The  effects  of  the  ionosphere  on  electromagenetic  waves  are  of  great 
interest  to  both  the  scientific  and  military  communities.  Radio  astronomical 
*  observations,  high  frequency  (HF)  radio  commimication,  and  over-the- 

horizon  (OTH)  radar  are  just  a  few  of  the  applications  upon  which  variations 
in  ionospheric  conditions  have  great  impact.  The  defining  characteristic  of 
ionospheric  behavior  is  the  density  of  free  electrons.  Electron  density 
fluctuates  widely  both  in  time  and  location  due  to  solar  activity,  geomagnetic 
activity  and  a  number  of  other  causes. 

The  current  method  of  measuring  electron  density,  which  has  been 
around  for  nearly  50  years,  is  to  transmit  a  frequency-swept  radio  signal  into 
the  atmosphere  and  record  the  reflection  at  the  ground  site.  While  this 
method  is  fast  and  accurate,  the  electron  density  profile  it  produces  is 
applicable  only  to  that  portion  of  the  ionosphere  immediately  above  the  site. 
In  order  to  produce  reliable  performance  in  HF  communication  or  OTH 
radar,  the  electron  density  must  be  known  over  global  areas.  The  20  or  so 
ionosonde  stations  currently  in  operation  cannot  provide  the  required 
\  geographic  coverage,  so  new  methods  are  being  explored.  One  possibility  for 

obtaining  truly  global  measurement  is  to  deploy  orbiting  ionosondes,  but  the 
power  and  mass  requirements  for  active  sensors  in  space  currently  make  this 
approach  impractical.  This  thesis  presents  a  possible  method  of  finding  the 
ionospheric  electron  density  over  a  large  area  using  passive  sensors. 


1 


The  sensors  in  this  case  are  the  Naval  Postgraduate  School's  (NFS) 

Middle  Ultraviolet  Spectrograph  (MUSTANG),  and  the  Naval  Research 

Laboratory's  (NRL)  High  Resolution  Airglow  and  Aurora  Spectrograph 

(HIRAAS).  The  method  is  the  measurement  and  analysis  of  the  ultraviolet 

(UV)  airglow  from  Earth's  upper  atmosphere.  Airglow  is  the  radiation  i 

emitted  from  atoms  and  molecules  which  are  in  excited  states  due  to 

photoexdtation,  electron  impact  excitation,  or  chemical  reaction.  HIRAAS 

was  designed  to  measure  the  emissions  of  atomic  oxygen  between  500  A  and 

15(X)  A,  and  will  not  be  discussed  in  this  thesis.  MUSTANG  was  designed  to 

observe  the  wavelength  region  from  1800  A  to  3400  A  in  which  the  UV 

emissions  from  molecular  nitrogen  (N2)  and  nitric  oxide  (NO)  are 

predominant.  The  purpose  in  studying  these  two  molecules  is  that  N2  is  the 

most  abundant  species  in  the  atmosphere  and  plays  an  important  part  in 

many  atmospheric  phenomena,  while  NO,  although  only  a  trace  constituent, 

is  a  strong  source  of  ions  in  the  lower  ionosphere.  Calculation  of  electron 

density  begins  with  the  measurement  of  emitted  intensities  which,  in 

principle,  allows  inference  of  the  volume  density  of  the  emitting  species. 

Once  the  densities  are  known,  models,  such  as  those  developed  by  Cleary 
(1986),  of  the  chemical  relationships  between  neutral  species  and  ions  can  be 
used  to  find  ion  and  electron  densities. 
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A.  THESIS  GOALS 

The  work  presented  in  this  thesis  has  two  goads:  to  establish  the 
sensitivity  calibration  for  the  MUSTANG  spectrograph  in  the  1800  A  to  2000 
A  wavelength  range,  and  to  determine  emission  profiles  for  N2  and  NO  by 
fitting  synthetic  spectra  to  observed  data. 

B.  THESIS  OUTLINE 

This  thesis  contains  five  chapters.  Chapter  n  presents  introductory 
information  on  topics  relevant  to  ionospheric  spectroscopy.  These  topics 
include  the  structure  of  Earth's  atmosphere  and  ionosphere,  molecular 
spectroscopy,  and  the  mechanisms  of  UV  emission.  Also  in  Chapter  n  is  a 
discussion  of  the  creation  and  utilization  of  synthetic  spectra. 

Chapter  in  describes  the  experiment  and  gives  a  brief  description  of  the 
MUSTANG  instrument.  The  vacuum  chamber  and  the  procedure  for 
vacuum  calibration  of  the  instrument  are  also  presented. 

Analysis  of  the  MUSTANG  data  and  results  of  fitting  with  synthetic 
spectra  are  given  in  Chapter  IV.  Each  of  25  observed  spectra  is  approximated 
with  a  synthetic  spectrum  in  the  wavelength  range  1850  A  to  2070  A.  The 
scaling  factor  which  produces  the  best  fit  for  each  spectrum  is  then  used  to 
determine  a  colunrn  emission  profile.  Conclusions  and  suggestions  for 
further  research  follow  in  Chapter  V. 
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IL  BACKGROUND 


The  purpose  of  this  chapter  is  to  provide  a  brief  introduction  to  the 
concepts  which  are  important  in  probing  the  ionosphere  with  a  rocket-borne 
spectrograph.  A  thorough  treatment  of  any  one  of  these  ideas  would  be  a 
thesis  in  itself,  so  the  material  presented  herein  is  necessarily  brief. 

A.  EARTH'S  ATMOSPHERE 

1.  The  Neutral  Atmosphere 

Any  desription  of  the  earth's  atmosphere  is  complicated  by  the 
fluid  nature  of  its  thermal,  chemical,  and  electrical  characteristics.  The  two 
most  common  characterizations  of  the  atmosphere  are  by  temperature  and  by 
chemical  composition.  Characterization  by  temperature  divides  the 
atmosphere  into  four  layers  as  shown  in  Figme  2-1.  The  first  layer,  the 
troposphere,  extends  from  Earth's  siuface  to  roughly  15  km  at  mid-latitudes. 
Convective  tiubvilence  is  the  major  physical  process  in  the  troposphere  and  is 
caused  by  the  heating  of  the  atmosphere  from  below  by  the  warm  surface. 
Tropospheric  temperahire  decreases  fairly  evenly  with  altitude  to  about  210  K 
at  the  tropopause. 

Above  the  tropopause  lies  the  stratosphere,  a  layer  about  25  km 
thick  at  mid-latitudes.  The  stratosphere  exhibits  a  positive  temperature 
gradient  due  to  absorption  of  solar  ultraviolet  radiation  near  the  base  of  the 
layer.  Ozone  is  a  particularly  strong  absorber  of  ultraviolet  radiation  in  the 
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2000  A  to  3000  A  wavelength  range.  The  temperature  incre2ise  stops  at  the 
stratopause  at  an  altitude  of  approximately  40  to  50  km. 

The  altitude  region  from  about  50  km  to  80  km  is  called  the 
mesosphere  and  is  another  region  of  negative  temperatiure  gradient. 
Temperature  in  the  mesosphere  is  controlled  largely  by  convection,  with  the 
warmer  stratosphere  below  affecting  the  layer  much  as  the  planet's  surface 
affects  the  troposphere.  The  atmosphere  reaches  its  lowest  kinetic 
temperature  <'»f  about  180  K  at  the  mesopaus?  near  80  km  altitude. 

Extending  from  the  mesopause  to  the  ill-defined  outer  edge  of  the 
atmosphere,  the  thermosphere  is  marked  by  a  strong  positive  temperature 
gradient  near  its  base.  There  are  several  factors  which  cause  the  gradient, 
notably  the  very  low  density  of  the  ga'-.es,  the  extremely  small  number  of 
triatomic  molecules,  and  heating  of  the  atmosphere  from  above  by  short 
wavelength  solar  radiation.  With  the  atmosphere  heated  from  above  rather 
than  from  below,  convection  ceases  and  kinetic  temperature  rises  to  a  level 
dependent  on  solar  activity.  The  high  temperatures  of  the  thermosphere  are 
the  result  of  the  lack  of  heat  removal  processes,  other  than  radiation,  at  low 
particle  density.  Beyond  400  km  to  500  km  altitude,  the  temperature  stops 
rising  and  becomes  independent  of  altitude. 

Characterization  by  chemical  composition  leads  to  a  three-layer 
description  of  the  atmosphere.  Temperatiues  and  densities  in  the 
troposphere,  stratosphere,  and  mesosphere  are  sufficiently  high  that 
convective  mixing  results  in  uniform  distribution  of  all  Cie  atmospheric 
gases.  This  region  of  uniform  distribution  is  called  the  homosphere  and  has  a 
composition  of  approximately  78%  molecular  nitrogen,  21%  molecular 
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oxygen,  and  1%  other  geises.  The  upper  limit  of  the  homosphere,  called  the 
turbopause,  is  marked  by  stable  temperature  gradients  and  the  emergence  of 
molecular  diffusion  rather  than  convective  mixing  as  the  primsuy  physical 
process.  Molecular  diffusion  and  the  subsidence  of  mixing  above  the 
turbopause  cause  the  atmospheric  species  to  stratify  according  to  mass, 
resulting  in  a  region  called  the  heterosphere.  In  the  heterosphere,  nitrogen 
remains  largely  m  molecular  form  while  photodissociation  of  molecular 
oxygen  produces  significant  amounts  of  atomic  oxygen  above  100  km. 

Above  approximately  500  km,  lies  the  exosphere  which  is  made 
mostly  of  helium  and  hydrogen  at  high  kinetic  temperature  and  which  forms 
an  indistinct  boundary  between  Earth’s  atmosphere  and  space  somewhere 
above  1000  km.  The  exosphere  is  marked  by  densities  so  low  that  particle 
trajectories  are  no  longer  the  random  paths  caused  by  collisions  with  other 
particles,  but  become  the  parabolic  trajectories  of  free-falling  bodies. 

2.  The  Ionosphere 

Although  most  of  the  atmosphere  is  made  of  neutral  atoms  and 
molecules,  portions  of  the  upper  atmosphere  have  ions  and  free  electrons 
present  in  sufficient  numbers  to  create  a  plasma  which  affects  radiowave 
propagation.  Such  regions  compose  the  ionosphere.  The  ionosphere  is  a 
weak  plasma  in  that  even  the  m^udmum  free  electron  concentration  is  three 
orders  of  magnitude  smaller  than  the  concentration  of  neutral  species.  The 
ionosphere  is  categorized  into  layers,  labelled  D,  E,  and  F,  according  to  altitude 
and  formation  process. 

The  production  of  charged  particles  in  the  upper  atmosphere  is  due 
primarily  to  ionization  of  neutral  gas  atoms  and  molecules  by  solar 
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ultraviolet  and  x-radiation.  The  creation  of  layers  of  electron  density  can  be 
explained  by  stud3dng  a  representative  reaction.  The  ionization  of  molecule 
Q  can  be  written  as 

Q  +  hv  — >  Q‘*‘  +  e'  . 

The  rate  of  formation  of  free  electrons  depends  on  the  product  of  the 
concentration  of  Q  and  the  flux  of  photons.  Particle  density  falls 
exponentially  while  the  flux  of  ionizing  photons  asymtotically  approaches  a 
maximum  as  altitude  increases,  resulting  in  peak  in  electron  production  at 
some  altitude  (see  Figure  2-2).  This  electron  production  peak  creates  an 
ionospheric  layer. 

The  lowest  layer  is  the  D  region  whirh  extends  from  about  50  km  to 
90  km.  The  primary  source  of  ions  in  the  D  region  is  photoionization  of  NO 
by  solar  Lyman-a  radiation,  written  as 

NO  +  hv  — >  NO"*"  +  e*  . 

Da5rtime  electron  densities  are  between  10^  and  10^  cm'^.  Relatively  high 
particle  densities  lead  to  rapid  recombination  of  ions  and  electrons  and  cause 
the  the  D  region  to  disappear  at  night.  Other  significant  sources  of  ions  are 
three-body  attachment  reactions,  e.g. 

O2  +  e'  +  M  — >  02"  +  M , 
zmd  radiative  attachment,  such  as 

O2  +  e’  — >  Q2'  +  hv  . 

The  major  loss  processes  in  the  D  region  are  collisional  detachment, 
photodetachment,  and  mutual  neutralization.  The  last  process,  exemplified 
by  the  reaction  +  B"  — >  A  +  B  ,  indicates  the  presence  of  large  niunbers 
of  negative  ions  in  the  layer. 
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Figure  2-2  Ionospheric  Layer  Formation 
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Between  90  km  and  140  km  lies  the  E  region  which  is  formed 
mostly  by  the  absorption  of  far-  and  extreme  ultraviolet  (FUV  and  EUV)  and 
soft  x-rays  by  diatomic  2md  monatomic  oxygen.  NO'*’  and  Q2'*'  are  the  most 
abundant  ions  in  the  E  region,  with  a  smaller  concentration  of  O'*'.  T)rpical 
electron  density  in  the  E  region  is  10^  cm*^.  The  principal  loss  mechanism  is 
dissociative  recombination,  for  example 

NO+  +  e-  ->  N  +  O  . 

Like  the  D  layer,  the  E  layer  tends  to  disappear  at  night,  except  at  high 
latitudes  where  energetic  particle  bombardment  of  the  upper  atmosphere 
causes  continued  ionization  after  sunset. 

The  O'*"  concentration  increases  with  altitude  until  it  becomes  the 
dominant  ion  in  the  F  region.  The  F  region  is  subdivided  into  FI  and  F2 
layers,  the  FI  layer  occup)dng  the  140  km  to  200  km  altitude  range,  and  the  F2 
layer  above  200  km.  The  presence  of  an  FI  layer  as  a  discemable  entity  is 
dependent  on  solar  activity.  When  it  is  present,  typical  electron  density  in  the 
FI  layer  is  2  x  10^  cm'^.  The  principal  ions  in  the  FI  layer  are  O'*'  and  NO'*'.  In 
the  absence  of  ionizing  radiation,  these  ions  quickly  disappear  via  charge 
exchange  and  dissociative  recombination.  Thus,  the  FI  layer  disappears  at 
night.  The  uppermost  layer  of  the  ionosphere,  the  F2  layer,  has  the  largest 
free  electron  density  in  the  ionosphere.  Daytime  density  is  of  the  order  of  2  x 
10^  cm"®  at  a  nominal  altitude  of  300  km.  Absorption  of  EUV  radiation  by 
atomic  oxygen  is  the  major  ionization  process  in  the  F2  layer,  and  O'**  the 
most  abundant  ion.  Part  of  the  reason  for  the  relatively  large  electron  density 
in  the  F2  layer  is  the  fact  that  ion-molecule  charge  exchange  replaces 
dissociative  recombination  as  the  major  loss  mechanism.  This  causes  the  loss 
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rate  of  ions  to  diminish  as  altitude  and  particle  density  decrease.  Relatively 
low  densities  and  bombardment  by  energetic  particles  allow  the  F2  region  to 
persist  at  night  with  electron  densities  about  five  percent  of  daylight  levels. 

B.  DIATOMIC  MOLECULAR  SPECTRA 

The  radiation  spectrum  of  a  molecule  is  a  measure  of  the  change  in  its 
energy  as  it  changes  state.  The  energy  state  of  a  diatomic  molecule  can  be 
described  by  four  quantities.  The  first  quantity  is  the  kinetic  energy  of  the 
molecule  due  to  the  velocity  of  its  center  of  mass.  This  energy  is  linearly 
proportional  to  the  temperature  of  the  molecule,  and  it  can  have  any  of  a 
continuiun  of  values.  When  the  temperature  of  the  upper  atmosphere  is 
said  to  be,  for  example,  1400  K,  it  is  equivalent  to  saying  the  average  kinetic 
energy  of  the  molecules  and  atoms  in  the  region  is 

1400  •  k  =  1.93  X 10'^^  J  =  0.121  eV  (2.1) 

where  k  =  1.381  x  10'^^  J/K  is  Boltzmann's  constant.  Energies  of  the  order  of 
0.1  eV  correspond  to  infrared  radiation,  which  is  not  examined  in  this  thesis. 
The  three  other  quantities  which  define  the  energy  state  of  a  diatomic 
molecule  are  related  to  the  internal  struchu'e  of  the  molecule,  sp>ecifically  its 
electronic,  vibrational,  and  rotational  states.  In  accordance  with  the  laws  of 
quantum  mechanics,  these  terms  can  have  only  discrete  values.  The 
quantization  of  energy  values  is  reflected  in  the  radiation  produced  when  the 
values  change.  Such  radiation  is  not  continuous,  but  forms  spectra  with 
emission  peaks  at  wavelengths  corresponding  to  the  energy  of  the  allowed 
transitions  for  the  radiating  molecule. 
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1.  Energy  Levels 

The  subject  of  molecular  spectrospcopy  is  treated  in  several  texts, 
and  the  following  explanation  is  taken  largely  from  that  of  Herzberg  (1950). 
The  energy  state  of  a  molecule  can  be  found  by  modelling  its  motions.  The 
simplest  model  of  a  diatomic  molecule  is  one  of  two  masses  connected  by  a 
massless  rigid  rod.  In  addition  to  translation  of  its  center  of  mass,  there  are 
two  possible  motions  for  the  molecule:  rotation  about  the  center  of  m2iss,  and 
vibration  along  the  intemuclear  axis. 

Solution  of  the  Schrodinger  equation  for  a  rigid  rotator  leads  to 
energy  levels  given  by 


EC)* 


(2.2) 


where  J  (=0,1,2,...)  is  the  rotational  quantum  number  and  h  is  Planck's 
coi\stant.  The  term  |ir^  is  the  moment  of  inertia  for  the  dumbbell  model. 
When  the  energy  of  a  state  is  measured  relative  to  the  groimd  state  it  is  called 
a  term  value.  The  rotational  term  value  is  the  rotational  energy  divided  by 
the  constant  he,  and  is  given  the  symbol  pq).  To  calculate  the  vibrational 
energy  levels,  the  rod  in  the  model  is  replaced  by  a  spring,  resulting  in  an 
harmonic  oscillator.  The  solution  of  Schrodinger's  equation  for  a  simple 
harmonic  oscillator  is  well  known  and  results  in  energy  levels  given  by 


(2^) 


where  v  is  the  vibrational  quantum  number,  p  the  reduced  mass,  and  k  the 
"spring  constant"  or  force  constant  for  the  interaction  between  the  atoms. 
Analogously  to  the  rotational  term  value,  the  vibrational  term  value,  G(v),  is 
E(v)  divided  by  he.  Since  real  molecules  are  neither  rigid  rotators  nor  true 
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harmonic  oscillators,  the  expressions  derived  above  for  term  values  are 
approximations.  Nonetheless,  the  expressions  are  conceptually 
straightforward  and  illustrate  the  origin  of  multiple  spectral  lines. 

There  is  a  third  way  in  which  the  energy  state  of  a  molecule  can 
change,  namely  by  making  a  transition  between  electronic  energy  levels.  As 
is  true  for  individual  atoms,  the  electrons  in  a  molecule  have  a  total  orbital 
angular  momentum,  denoted  L.  In  general,  L  will  precess  about  the 
intemuclear  axis,  but  the  component  of  L  along  the  axis  is  a  constant  of  the 
electronic  motion,  and  is,  therefore,  a  useful  measure  of  the  electronic  state  of 
the  molecule.  This  axial  component  of  L,  called  Ml,  is  quantized  and  can 

have  values  L,  L-1,  L-2,...,-L.  The  magnitude  of  Ml  is  given  the  s5mribol  A,  so 
that  A  =  |Ml|  =  0,1,2,. ..,L.  The  values  of  A  for  molecular  states  are  given  Greek 

letters  corresponding  to  the  Roman  letters  used  to  designate  orbital  angular 
momentum  for  atomic  states.  States  £,  11,  and  A  correspond  to  A  =  0,l,and  2. 
In  addition  to  their  orbital  angular  momentum,  the  electrons  in  a  molecule 
have  spin  angular  momentum  which  generates  a  resultant  spin,  S.  The  spin 
will  precess  around  the  intemuclear  axis  in  the  presence  of  a  magnetic  field, 
but,  as  with  the  orbital  angular  momentum,  the  component  of  S  along  the 
axis  is  a  constant  of  the  electron  motion.  This  component,  labelled  £,  can  be 
oriented  either  parallel  or  anti-parallel  to  the  intemuclear  axis,  so  there  are 
two  values  of  the  spin  quantum  number  in  the  presence  of  a  magnetic  field. 
When  A  s  0,  there  is  no  internal  magnetic  field  in  the  molecule,  and  there  is 
only  one  value  for  the  axial  component  of  spin.  When  A  is  not  0,  the  orbiting 
electrons  create  an  internal  magnetic  held,  giving  rise  to  2S  values  of  spin. 

All  together  then  there  are  2S  +  1  possible  values  of  spin  angular  momentum 
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for  a  molecule.  The  quantity  (2S  +  1)  is  called  the  multiplicity  of  the  state  of 
the  molecule.  The  term  value  for  electronic  energy  is  given  the  symbol  T®. 
The  values  of  Te  are  determined  by  laboratory  measurement  and  have  been 
compiled  for  most  diatomic  molecules  in  Huber  and  Herzberg  (1979). 

In  spectroscopic  notation,  the  multiplicity  of  a  state  is  written  as  a 
left  superscript  to  the  orbital  angular  momentum  sjmibol.  For  example,  a 
state  having  spin  quantum  number  0  (£  +  1  =  1)  and  orbital  angular 
momentum  quantum  number  A  =  1  is  written  ^11.  The  electronic 
configuration  of  the  molecule  is  also  given  a  letter  designation.  X  is  used  to 
identify  the  grotmd  state,  with  succeeding  excited  states  given  S3anbols  A,  B,  C, 
or  a,  b,  c  depending  on  whether  the  excited  state  has  multiplicity  the  same  as 
(upper  case)  or  different  from  (lower  case)  the  ground  state.  A  complete 
description  of  molecular  spectroscopy  notation  may  be  foimd  in  Herzberg 
(1950). 

2.  Spectral  Lines 

The  relationships  between  electronic,  vibrational,  and  rotational 
energy  levels  are  shown  in  Figure  2-3.  Except  in  the  case  of  a  single  atom  with 
Ss:0,  the  radiation  h'om  a  change  in  energy  is  not  a  single  line,  but  a  multiplet 
of  lines  which  are  closely  spaced  in  wavelength.  In  diatomic  molecules,  the 
many  possible  combinations  of  rotational  and  vibrational  states  lead  to  multi- 
line  spectra  even  for  transitions  between  electronic  states  with  multiplicities 
of  one.  Figure  2-4  shows  the  effect  of  the  vibrational  and  rotational  terms  on 
the  spectrum  of  the  electronic  transition  from  the  first  excited  state  to  the 
ground  state  of  nitric  oxide.  Note  in  the  hgure  that  the  spacing  between  lines 
can  be  fractions  of  an  Angstrom,  corresponding  to  energy  ditferences  of  the 
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Figure  2-3  Relationship  Between  Electronic,  Vibrational,  and 

Rotational  Energy  Levels  (Eisberg  and  Resnick,  1985) 
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High  Resolution  Spectrum  Showing  Vibrational  and 
Rotational  Structure  (Qeary,  private  communication) 
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order  of  lO"^  eV  between  rotational  energy  levels.  When  measured  with  an 
instrument  such  as  MUSTANG,  which  has  a  wavelength  resolution  of 
approximately  10  A,  the  many  lines  appear  as  one  broad  peak.  Compare 
Figure  4-2. 


The  wave  number,  o,  of  a  spectral  line  is  determined  by  the 
difference  in  energies  between  the  initial  and  hnal  states,  and  is  given  by 


Ei-Ef 

^=“h^  • 

From  the  preceeding  discussion  on  term  values,  the  term  value,  T,  for  a 
given  state  may  be  written 


(2.4) 


T  =  Te  +  G(v)  +  Fa)  •  (2.5) 

Using  the  convention  of  denoting  the  upper  state  of  a  transition  with  a  single 
prime  and  the  lower  state  with  a  double  prime,  the  wave  number  for  a  given 
transition  can  be  written 


o  =  T'  -  T'  =  T’e  -  T'e  +  G(v’)  -  G(v")  +  FQ’)  -  F(J")  .  (2.6) 

The  group  of  transitions  for  which  the  upper  and  lower  electronic  energies 
are  fixed  is  referred  to  as  a  molecular  bai  d  system.  This  thesis  is  concerned 
with  the  LBH  bands  of  molecular  nitrogen  and  the  y,  5,  and  E  bands  of  nitric 
oxide.  For  each  of  these  band  systems,  the  lower  electronic  state  is  the  groimd 
state. 
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C  TRANSITION  MECHANISMS  FOR  N2  AND  NO 


1.  Selection  Rules 

The  airglow  spectrum  in  the  1800  A  to  2000  A  wavelength  range  is 
composed  of  the  LBH  bands  of  N2  and  the  y.  5,  and  e  bands  of  NO.  The 
specific  transitions  involved  are: 


y: 

A^I'^ 

<->  X^Hr 

6: 

c^n 

<->  X^Hr 

e : 

<->  X^Hr 

in  NO  and 

a^Hg  <->  X^I+g 


in  molecular  nitrogen  for  LBH.  The  y,  5,  and  e  bands  of  NO  arise  from 
allowed  transitions,  but  the  LBH  transition  is  forbidden.  The  term 
'forbidden'  does  not  indicate  that  a  particular  transition  never  occurs. 

Instead,  it  means  that  the  selection  rules  governing  the  occurrence  of  such  a 
transition  via  electric  dipole  radiation  are  violated.  For  a  brief  explanation  of 
selection  rules,  consider  the  case  of  a  one-electron  atom  with  a  dipole 
moment  p  *  -  er,  where  e  is  the  electron  charge  and  r  is  the  vector  from  the 
center  of  the  nucleus  to  the  electron.  The  probability  per  unit  time  that  the 
atom  will  emit  a  photon,  the  transition  rate,  can  be  written 


R« 


IbTC^P^fi 


/ 


where  v  is  the  frequency  of  the  photon.  The  term  pf^  is  called  the  matrix 
element  of  the  electric  dipole  moment  between  the  initial  and  final  states, 
and  is  defined  by  the  equation 
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|\l/*fer\jfidT 


(2.7) 


Pfi  = 


where  and  f  are  the  eigenfunctions  describing  the  initial  and  final  energy 
levels  and  dt  indicates  integration  over  all  space.  The  selection  rules  are  a  set 
of  conditions  on  the  quantiun  numbers  associated  with  yi  and  yf,  which  lead 
to  a  value  of  zero  for  when  those  conditions  are  violated  (Eisberg  and 

Resnick,  1985).  Violation  of  the  selection  rules  leads  to  a  transition  rate  of 
zero,  hence,  a  "forbidden"  transition.  The  rules  for  multi-electron  atoms  and 
molecules  are  derived  in  a  similar  maimer,  but  with  vastly  more  complicated 
wavefunctions. 

The  designation  of  the  LBH  band  trcursition  as  forbidden  means 
that  it  is  not  driven  by  a  change  in  the  electric  dipole  moment  of  the  nitrogen 
molecule.  In  particular,  the  LBH  transition  violates  a  rule  which  requires  the 
change  in  the  sum  of  the  orbital  angular  ixiomenta  of  the  electrons  in  the  N2 
molecule  to  be  AL  =  ±  1  for  electric  dipole  radiation.  There  is,  however, 
another  rule  which  says  that  the  electron  angular  momentum  may  change  by 
AL  ±  2  and  produce  radiation  via  a  change  in  the  electric  quadrupole  or 
magnetic  dipole  moments  of  the  molecule.  It  is  this  rule  which  leads  to  LBH 
radiation.  Because  of  these  selection  rules,  the  excitation  mechanism  for  LBH 
is  not  photoexcitation  but  energetic  photoelectron  impact.  In  contrast,  the 
nitric  oxide  bands,  which  come  from  allowed  transitions,  are  driven  by 
photoexcitation  and  emit  via  electric  dipole  radiation. 

2.  Volume  Emission  Rates 

The  population  of  a  state  depends  on  the  electronic,  vibrational, 
and  rotatioiuil  term  values  which  define  it.  An  accurate  expression  for  the 
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population  requires  that  corrections  be  made  to  the  expressions  for  term 
values  derived  from  the  simple  model  of  the  diatomic  molecule. 

The  first  correction  accoimts  for  the  anharmonidty  of  real 
vibrat'Dnal  motion.  The  corrected  vibrational  term  value  is  written 

G(v)  =  (De(v  + -  a>ex^.(v  + .  (2.8) 

(i>e  and  Xe  are  constcmts  of  the  molecule  and  are  listed  in  Huber  and  Herzberg 
(1979).  Corrections  to  the  rotational  term  values  involve  the  multiplicity  of 
the  state  and  are  more  complicated  than  those  for  the  vibrational  terms.  For 
singlet  states,  the  corrected  rotational  term  value  is 


Fa)  =  Bvjg  + 1)  -  Dvj^a  + 1)^  . 

(2.9) 

The  Bv  and  Dy  terms  depend  on  the  vibrational  state  of  the  molecule  and  are 

found  from  the  expressions 

By  *  Be  -  ae(v  +  |)  +  higher  order  terms. 

and 

(2.10) 

Dy  =  De  -  Pe(v  +  |)  +  higher  order  terms 

(2.11) 

The  constants  Be,  De,  o^,  and  again  found  in  Huber  and  Herzberg. 
When  the  multiplicity  of  the  state  is  greater  than  one,  the  rotational  term 
value  depends  strongly  upon  the  coupling  between  orbital  and  spin  angular 
momenta.  Hund  (1927)  classified  five  coupling  schemes  based  on  the 
strength  of  the  coupling.  These  schemes  will  not  be  explored  in  this  thesis, 
but  they  are  discussed  in  Cleary  (1985)  and  Bosserman  (1989). 
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The  volume  emission  rate  of  a  specific  spectral  line  is  the  product 
of  the  population  of  the  upper  state  to  which  the  line  corresponds  and  the 
coefficient  for  spontaneous  emission  between  the  upper  and  lower  states. 
Using  the  notation  VE  for  volume  emission  rate,  N  for  the  fraction  of 
molecules  in  the  upper  state,  A  for  the  spontaneous  emission  coefficient,  and 
n  to  denote  electronic  state,  the  emission  rate  can  be  written 


E(n’v'N7'p',  n"v"N"J"p")  =  N(n’v'N'J'p')  A(nV'N'J’p',  n"v"N"J"p")  ,  (2.12) 


where  n,v,N,J,  and  p  describe  the  electronic,  vibrational,  and  rotational  state 
according  to  Himd's  models.  For  an  upper  level  defined  by  n’  and  v’,  the 
fractional  population  of  the  level,  in  steady  state,  is  found  by  dividing  the 
production  rate,  P(n'v'),  by  the  loss  rate  for  that  level.  In  both  the  LBH  and 
the  nitric  oxide  bands,  the  loss  rate  for  the  upper  level  is  just  the  sum  of  the 
spontaneous  emission  rates  to  all  possible  lower  vibrational  levels. 
Quenching  can  be  ignored  for  LBH  above  95  km  (Vallance  Jones,  1974).  The 
fractional  population  is  then  written 


N(nV) 


P(n’v') 

X^v’v’ 

v" 


(2.13) 


For  the  nitric  oxide  bands,  the  production  rate  depends  on  the  solar  flux,  the 
absorption  cross-section  for  transition  to  the  upper  state,  and  on  molecular 
number  density.  Using  the  notation  of  Barth  (1965),  this  rate  can  be  written  as 

P(nV)  =  nF  ^  £.„  qo,.  ai« 

where  tcF  is  the  solar  flux,  X,  is  the  wavelength  of  the  incoming  photon,  f^^  is 


the  oscillator  strength  for  the  electronic  transition  from  groimd  state  X  to 
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upper  state  n,  and  qQ^.  is  the  Franck-Condon  factor  for  the  v=0  level  of  the 
ground  state  to  the  v'  level  of  the  excited  state.  Franck-Condon  factors  are 
explained  below.  In  the  case  of  the  LBH  bands,  the  excitation  rate  depends  on 
the  number  density  and  on  the  photoelectron  flux  and  electron  impact  cross- 
section  of  the  molecides.  The  production  rate  in  this  case  can  be  expressed  as 
follows: 

P(nV)  =  qov'Nj0(E)o(E)dE  ,  (2.15) 

where  N  is  the  number  density,  <t>(E)  denotes  the  flux  of  photoelectrons,  a(E)  is 
the  cross-section  for  an  X-to-n'  electronic  transition,  and  dE  indicates 
integration  over  the  range  of  photoelectron  energies.  The  difference  in 
excitation  processes  for  the  two  molecules  makes  calculation  of  the  excitation 
rate  for  NO  much  easier  than  that  for  N2. 

The  flux  of  photons  with  the  correct  energies  to  excite  the  NO 
bands  exhibits  little  variability  over  the  solar  cycle.  The  absorption  cross- 
sections  of  NO  for  these  photons  are  constants  and  have  been  measured.  On 
the  other  hand,  the  flux  of  the  photoelectrons  which  produce  the  upper  state 
of  the  LBH  transition  is  highly  variable.  The  primary  source  of 
photoelectrons  is  the  ionization  of  N2,  O2,  and  O  by  solar  EUV  radiation, 
particularly  the  He  304  A  emission.  Variations  in  the  sun's  output  become 
larger  as  the  wavelength  of  the  radiation  gets  shorter,  so  the  EUV  emissions 
which  generate  photoelectrons  fluctuate  more  than  do  the  MUV  emissions 
sustaining  NO  excitation.  Also,  these  photoelectrons  have  sufficient  kinetic 
energy  to  cause  multiple  ionizations.  The  secondary  photoelectrons  penetrate 
into  the  atmosphere  causing  more  ionization  and  losing  energy  in  collisions 
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until  they  are  thermalized.  The  energy  spread  of  this  cascade  of 
photoelectrons  must  be  modelled  by  some  form  of  random  distribution. 
Since  the  collision  cross-section  of  nitrogen  is  a  function  of  energy,  accurate 
calculation  of  the  LBH  excitation  rate  requires  knowledge  of  the  energy 
distribution  of  the  photoelectrons. 

When  the  excitation  rate  has  been  established,  the  population  of 
specific  states  must  be  determined.  To  make  this  determination,  the 
branching  ratio  into  the  available  vibrational  and  rotationcd  levels  must  be 
found.  For  given  electronic  levels,  the  branching  ratio  into  the  available 
vibrational  levels  is  fovmd  in  the  Franck-Condon  factors.  The  Franck- 
Condon  factors  represent  the  transition  probabilities  between  upper  and 
lower  vibrational  states.  When  the  wavefunctions  for  the  vibrational  states 
are  known,  the  Franck-Condon  factor  can  be  calculated  from 


qv’v"  = 


J\|/uer\|/jdx 


(2.16) 


where  v'  is  the  upper  vibrational  state  with  wavefunction  Vu/  «»nd  v"  is  the 
lower  vibrational  state  with  wavefunction  Yj.  When  the  wavefimctions  are 

not  known,  the  factors  for  a  band  progression  must  be  measured.  The  sum  of 
cdl  Franck-Condon  factors  for  a  fixed  value  of  v'  or  a  fixed  value  of  v"  is 
unity,  signifying  that  the  vibrational  transition  that  actually  occurs  must  be 
one  of  those  available.  When  the  Franck-Condon  factor  is  relatively  large, 
the  radiation  from  the  corresponding  vibrational  transition  is  strong.  This 
fact  is  important  in  the  fitting  of  synthetic  spectra  to  the  measured  data 
presented  in  Chapter  IV.  The  branching  ratio  from  a  particular  vibrational 
level  into  each  of  the  available  rotational  states,  J,  is  determined  by  the  Honl- 
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London  factors.  The  sum  of  the  Honl-London  factors  for  a  fixed  value  of  J’  is 
(2J'  +  1),  so  the  Honl-London  factor  for  a  specific  rotational  transition  must  be 
divided  by  (2J'  +  1)  to  give  a  probability. 

Assembling  Equations  (2.12),  (2.13),  and  (2.14)  or  (2.15)  as 
appropriate,  gives  a  Hnal  expression  for  the  volume  emission  rate  as 


E  .  N(nVNTp’)  A,.,.  SOT) 

E  -  r(nv)  N(n'v’)  2J' +  1 

V" 


(2.17) 


In  measuring  emission  rates,  accoimt  must  be  taken  of  the 
phenomenon  of  self-absorption.  Self-absorption  occurs  when  a  photon 
emitted  from  an  exdted  molecule  is  absorbed  by  a  molecule  of  the  same 
species.  Since  the  molecules  have  identical  energy  state  structures,  the  energy 
of  the  emitted  photon  will  match  exactly  one  of  the  transition  energies  in  the 
molecule  it  strikes,  and  the  photon  may  be  absorbed.  For  this  to  be  true,  the 
absorbing  molecule  must  be  in  the  same  state  as  was  the  emitting  molecule 
after  it  emitted  the  photon.  Since  the  vast  majority  of  molecules  in  the 
atmosphere  are  in  the  lowest  vibrational  state,  only  those  photons  from 
transitions  to  v"sO  are  likely  to  be  absorbed.  Note  that  since  the  LBH 
transition  is  driven  by  photoelectron  impact  rather  than  photoexcitation,  its 
emissions  are  not  subject  to  self-absorption. 


D.  SYNTHETIC  SPECTRA 

In  order  to  analyze  the  observed  spectra,  models  are  needed  which  create 
S3mthetic  spectra  based  on  theoretical  considerations  such  as  those  which  lead 
to  Equation  (2.17).  Such  models  have  been  developed  by  Cleary  (1986)  for  the 
y,  6,  and  e  bands  of  nitric  oxide  and  for  the  LBH  bands  of  N2.  The  basis  of  the 
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model  is  the  calculation  of  a  theoretical  volume  emission  rate  for  each 
rotational  line  in  a  molecular  band.  The  spectral  region  of  concern  is  divided 
into  wavelength  bins,  and  the  emission  rates  are  placed  in  the  appropriate 
bins  to  create  a  synthetic  spectrum. 

Before  a  spectrum  from  one  of  these  models  can  be  vised  for  comparison 
with  observed  data,  it  must  be  conditioned  to  reflect  the  characteristics  of  the 
measuring  instrument.  To  do  so,  the  S3mthetic  spectrum  is  convolved  with 
the  instrument  slit  function.  The  slit  function  is  the  frequency  domain 
transfer  fimction  of  the  spectrograph,  and  is  determined  by  measuring  the 
instrument  response  to  a  wavelength  delta  function.  For  the  MUSTANG  slit 
function  determination,  the  delta  functions  were  approximated  by  several  of 
the  spectral  lines  of  atomic  mercviry.  The  wavelength  resolution  of  the 
MUSTANG  was  found  to  be  10.6  A,  and  the  effect  of  convolution  on  the 
S5mthetic  spectra  can  be  seen  in  Figure  2-5.  As  noted  earlier,  the  resolution  of 
the  MUSTANG  is  too  coarse  to  show  the  rotational  structure  of  a  band. 

One  purpose  of  fitting  the  observed  spectra  with  synthetic  ones  is  to  find 
coluirm  densities  of  NO  at  each  altitude.  As  explained  previously,  the 
excitation  rate  for  NO  is  calculable,  so  the  measured  emission  rate  is  directly 
related  to  the  density  of  nitric  oxide.  The  relationship  between  the  synthetic 
fit  and  colvimn  density  is  straightforward.  The  intensity  recorded  by  the 
spectrograph  is  a  volume  emission  rate  integrated  over  the  instrument  line 
of  sight,  while  the  synthetic  spectrum  represents  the  volume  emission  rate 
from  a  single  molecule.  The  scale  factor  between  the  synthetic  and  observed 
spectra  is,  therefore,  equal  to  the  column  density  of  the  emitting  molecule. 
The  molecular  densities  then  become  inputs  to  other  models  which  lead  to 
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the  ultimate  goal  of  deducing  ionospheric  electron  densities.  The  use  of 
synthetic  spectra  to  infer  molecular  nitrogen  density  from  LBH  emissions  is 
not  as  direct  as  it  is  for  nitric  oxide.  While  the  scale  factor  for  the  LBH  spectra 
still  represents  column  emission  rate,  as  it  does  for  the  NO  spectra,  the  rate 
cannot  be  directly  converted  to  a  density  for  nitrogen.  Many  of  the 
parameters  that  affect  the  LBH  emission  rate  are  not  directly  measureable,  and 
the  relative  variations  among  them  are  unknown.  For  example,  the  density 
of  nitrogen  is  one  of  the  important  factors  in  determining  the  photoelectron 
flux  driving  the  emissions  by  which  we  seek  to  measure  nitrogen  density. 

The  coupling  between  variables  makes  it  necessary  to  assume  values  for 
certain  parameters,  and  to  use  those  values  to  generate  coherent  results. 
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Figure  2-5  Effect  on  Synthetic  Spectrum  of  Convolution  with  MUSTANG 
Slit  Function 
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ffl.  THE  EXPERIMENT 


A.  INTRODUCTION 

The  MUSTANG  mstrument  constituted  half  of  an  ionospheric  research 
package  that  was  launched  at  1700  GMT  on  30  March,  1990  from  White  Sands 
Missile  Range  aboard  a  NASA  Terrier  boosted  Black  Brant  sounding  rocket. 
The  other  half  of  the  payload  was  the  NRL  HIRAAS,  an  instrument  designed 
to  measure  UV  emissions,  primarily  from  O'*',  in  the  500  A  to  1500  A 
wavelength  region.  Latitude  and  longitude  of  the  launch  site  are  32.4°  N  and 
106.34°  W,  respectively.  At  the  time  of  launch,  the  solar  zenith  angle  was 
62.9°,  and  F10.7  was  186.9.  Total  flight  time  was  approximately  14  minutes, 
with  the  payload  reaching  a  maximum  altitude  of  320  kilometers.  An  active 
attitude  control  system  kept  the  payload  properly  oriented  for  data  collection 
throughout  the  flight.  For  MUSTANG,  the  correct  orientation  was  a  90° 
observation-zenith  angle  (horizontal  view  in  the  rocket’s  frame  of  reference) 
and  90°  solar-observation  angle.  This  orientation  results  in  a  view  through 
the  limb  of  the  earth's  atmosphere.  Such  a  view  was  necessary  both  to 
prevent  saturation  of  the  detector  by  the  very  bright  airglow  in  the  nadir 
direction,  and  to  enable  detection  of  the  NO  and  N2  emissions  between  18(X) 

A  and  3000  A  which  are  almost  completely  absorbed  by  O2  O3  in  the  nadir 
view.  The  MUSTANG  instrument  has  no  inherent  recording  capability,  so 
the  data  were  telemetered  to  a  recording  station  on  the  ground.  The  payload 
section  of  the  rocket  was  recovered  after  a  parachute-retarded  final  descent. 
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B.  THE  INSTRUMENT 

The  MUSTANG  mstrument  is  a  space-quaMed  1/8  meter  Ebert-Fastie 
spectrograph  with  an  off-axis  telescope.  A  photodetector  consisting  of  an 
image  intensifier  opticzdly  coupled  to  a  linear  photodiode  array  is  positioned 
at  the  exit  focal  plane  of  the  spectrograph.  The  voltage  output  from  each  of 
the  512  photodiodes  in  the  array  is  proportional  to  the  intensity  of  radiation  at 
the  wavelength  to  which  the  diode  corresponds.  The  limits  of  wavelength 
coverage  are  hrom  1800  A  to  3400  A  which  is  the  middle  ultraviolet  portion  of 
the  electromagnetic  spectrum.  The  spectral  resolution  of  MUSTANG  is  10.6 
A;  and,  to  avoid  anomalovis  data  spikes,  the  radiation  at  a  given  wavelength 
was  sampled  by  three  adjacent  pixels,  each  of  which  has  a  coverage  of  3.133  A. 
For  a  complete  description  of  the  detector,  see  Clayton  (1990)  or  Andersen 
(1990). 

C  VACUUM  CALIBRATION 

Calibration  of  the  spectrograph  in  the  wavelength  range  2000  A  to  3400 
A  was  accomplished  by  Andersen  (1990)  in  his  analysis  of  the  OH  2470  A  and 
OI  2972  A  features.  An  extension  of  that  calibration  to  wavelengths  below 
2000  A  was  done  to  enable  analysis  of  the  LBH  bands  of  N2  from  1850  A  to 
2070  A.  Such  an  extension  required  the  removal  of  oxygen  from  the  path 
between  the  UV  calibration  source  and  the  spectrograph  in  order  to  preclude 
absorption  of  the  signal  by  the  O2  Shumann-Runge  bands.  Oxygen  absorption 
was  eliminated  by  performing  the  calibration  in  a  vacuum  chamber  which 
had  been  evacuated  and  refilled  with  nitrogen. 
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1.  Vacuum  Chamber 

The  vacuum  chamber,  shown  schematically  in  Figure  3-1,  was  built 
from  a  1.5  m  long  section  of  PVC  pipe  12  inches  in  diameter.  One  end  of  the 
chamber  was  sealed  with  an  end  cap  through  which  ports  were  drilled  and 
tapped  for  the  vacuum  and  nitrogen  backfill  lines.  A  vacuum  gauge  was 
fitted  to  the  nitrogen  line.  The  other  end  of  the  pipe  was  fitted  with  a  flange 
assembly  containing  an  O-ring  groove  and  an  attaching  ring  to  which  a  12.5 
mm-thick  aluminum  plate  was  bolted.  A  hole  in  the  center  of  the  aluminum 
plate  was  covered  by  an  110  mm  diameter  feed-through  plate  through  which 
a  37-pin  ribbon  cable  and  three  single  wires  were  passed.  The  ribbon  cable 
provided  coimection  to  the  spectrograph,  and  the  wires  provided  power  to  a 
deuterium  discharge  lamp.  The  deuterium  lamp  is  an  ultraviolet  spectral 
irradiance  standard  for  calibration  of  spectroradiometer  systems. 

Inside  the  chamber,  the  instrument  and  the  lamp  were  mounted 
on  a  4.75  mm  thick  aluminvim  table  measuring  1.4m  x  .285m.  The  lamp  was 
mounted,  at  a  height  corresponding  to  the  center  of  the  entrance  slit  in  the 
spectrograph,  on  a  micrometer-driven  translation  stage.  The  translation  stage 
allowed  the  precise  horizontal  positioning  of  the  lamp  required  to  get  a 
useable  signal  through  the  140-nm  width  of  the  slit.  The  distance  firom  the 
lamp  to  the  detector  was  1.83  m  and,  when  properly  positioned,  the  lamp 
output  was  sufficiently  strong  to  saturate  the  detector.  A  piece  of  black  card 
stock  with  a  pinhole  in  the  center  was  placed  over  the  end  of  the  spectrograph 
telescope  reducing  the  signal  intensity  to  a  non-saturating  level.  The  close 
proximity  of  the  lamp  to  the  wall  of  the  chamber  caused  concern  about  the 
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possibility  of  heat  and/or  UV  damage  to  the  chamber,  so  the  upper  section  of 
the  chamber  wall  in  the  vicinity  of  the  lamp  was  covered  with  aluminum 
foil. 

2.  Experimental  Procedure 

Spectra  were  recorded  under  several  conditions  once  the  vacuum 
chamber  apparatus  was  assembled.  In  all  cases,  the  deuterium  lamp  wais 
switched  on  and  allowed  ten  minutes  to  warm  up  before  any  spectra  were 
recorded.  Despite  the  warm  up  period,  a  slow  decrease  in  output  was  seen 
consistently  during  each  set  of  recordings.  Andersen  (1990)  attributes  this 
decrease  to  an  exponentially  decreasing  spectrograph  response  with  a  time 
constant  of  224  seconds.  The  decrease  in  mstrument  response  introduces  an 
uncertainty  of  ±8%  in  the  calibration. 

The  first  spectra  were  tsdcen  with  the  cover  plate  removed  from  the 
chamber  and  the  apparatus  open  to  the  laboratory  atmosphere.  Fifteen 
recordings  were  made  at  atmosphere,  each  consisting  of  fifty  individual 
spectra  averaged  together.  The  chamber  was  then  sealed  and  evacuated,  with 
cm  oil-sealed  mechanical  vacuum  pump,  to  a  pressure  of  approximately  three 
torr.  Power  was  removed  from  the  lamp  and  the  spectrograph  during 
pumping,  and  no  spectra  were  taken  imder  vacuum  conditions  in  order  to 
avoid  low-pressure  arcing  of  high-voltage  electronic  components  in  the 
spectrograph. 

Once  the  chamber  was  evacuated,  dry,  gaseous  nitrogen  was 
introduced  into  the  chamber  until  it  reached  atmospheric  pressure.  Power 
was  reapplied  to  the  lamp  and  spectrograph  and,  after  the  ten-minute  warm 
up,  fifteen  recordings  were  made.  This  cycle  was  repeated  five  times  for  a 
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total  of  six  purge/fill  cycles.  Since  one  atmosphere  exerts  approximately  760 
torr  pressure,  each  evacuation  of  the  chamber  reduced  the  concentration  of 
oxygen  to  3/760  ths  of  atmospheric  concentration.  The  amount  of  oxygen 
remaining  in  the  chamber  after  the  last  cycle  was  estimated  to  be 

(5|q)^  =  4x10-^5  (31) 

times  the  oxygen  concentration  in  the  laboratory  air. 

While  these  preliminary  calibration  data  were  being  studied,  the 
apparatus  was  allowed  to  remain  in  the  nitrogen-filled  chamber  at 
atmospheric  pressure  for  ten  days.  The  chamber  was  then  evacuated  and 
refilled  with  room  air.  Three  recordings  were  made  after  the  ten-minute 
lamp  warm  up  and  three  more  made  after  an  additional  20  minutes  of  lamp 
operation.  When  the  recording  was  complete,  the  chamber  was  opened,  and 
the  apparatus  removed.  After  a  week  in  the  laboratory  atmosphere,  six  more 
recordings  were  taken  to  provide  a  comparison  of  performance  before  and 
after  operation  in  the  vacuum  chamber. 

3.  Results 

Figure  3-2  clearly  shows  the  increase  in  detector  response  below 
2000  A  when  oxygen  is  removed  from  the  path.  The  O2  absorption  bands  at 
1972  A,  1947  A,  1924  A,  1902  A,  and  1882  A  can  be  seen  to  decrease  with 
successive  evacuations  of  the  chamber.  Although  the  system  was  purged  six 
times,  the  detector  response  did  not  increase  noticeably  after  the  third 
purge /fill  cycle,  at  which  point  the  chamber  oxygen  concentration  was 
decreased  by  more  than  seven  orders  of  magnitude.  The  system  response 
after  the  third  purge  was  used  as  the  short-wavelength  calibration  standard. 
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Figtire  3-2  Detector  Response  After  Oxygen  Purge 
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In  order  to  get  meaningful  numbers  from  the  data,  the  system 
ouput  must  be  converted  to  units  which  describe  atmospheric  radiation.  The 
unit  used  to  describe  radiation  is  the  Rayleigh,  which  is  defined  as  10^ 
photons/cm^^sec.  The  S5ratem  output  is  a  voltage  which  is  measured  in 
decimal  luiits  (du)  per  pixel.  Dividing  this  output  by  the  absolute  intensity,  in 
Rayleighs  per  Angstrom  (R/A),  of  the  deuterium  irradiance  standard  lamp 
pelds  a  calibration  parameter  with  units  of  (du/pixel)  per  (R/A).  When  the 
observed  MUSTANG  data,  which  have  units  of  (du/ pixel),  are  divided  by  the 
calibration  parameter,  the  result  is  spectra  with  units  of  (R/A).  These  spectra 
can  then  be  integrated  over  wavelength  to  find  an  emitted  intensity  in 
Rayleighs. 

The  calibration  performed  in  the  vacuum  chamber  differed  in  the 
physical  arrangement  of  the  equipment  from  that  done  by  Andersen(1990)  in 
atmosphere.  In  the  chamber,  the  spectrograph  was  placed  to  view  the  lamp 
directly  from  a  distance  of  1.83  m,  whereas  Andersen  positioned  the 
spectrograph  to  view  a  reflecting  screen  illuminated  by  the  lamp  from  a 
distance  of  2.00  m.  In  addition  to  the  different  physical  arrangements,  two 
important  equipment  considerations  led  to  a  decision  to  use  the  vacuum 
chamber  results  to  extend  Andersen's  calibration  rather  than  to  create  an 
entirely  new  calibration. 

The  first  consideration  is  the  fact  that  the  deuterium  lamp  output 
changes  with  lamp  age.  The  absolute  intensity  of  the  lamp  output  falls  as  the 
lamp  ages,  but  does  so  fairly  uniformly  across  the  wavelength  range  of 
interest  making  an  absolute  calibration  difficult,  but  having  little  effect  on  a 
relative  calibration.  The  second  consideration  arises  from  the  choice  of  PVC 
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as  a  vacuum  chamber  material.  PVC  pipe  was  chosen  for  the  chamber 
because  it  was  available  quickly  and  inexpensively^  and  the  chamber  was  not 
intended  for  use  at  high  vacuum.  Unfortimately,  three  da3rs  of  curing  and  80 
hours  of  continuous  purging  still  did  not  remove  from  the  chamber  all  the 
volatile  components  of  the  adhesives  used  to  seal  the  ends  of  the  pipe.  When 
the  apparatus  was  allowed  to  stand  in  the  chamber  for  ten  days,  the  volatile 
components  are  believed  to  have  formed  a  film  on  the  equipment.  The  effect 
of  the  film  can  be  seen  in  the  curve  labelled  'After  Vacuum  Chamber,  Film 
Present'  in  Figure  3-3.  After  a  week  in  the  open  laboratory,  the  film  had 
evaporated  and  detector  response  at  short  wavelengths  (pixels  0  to  92) 
resembled  pre-vacuum-chamber  response.  Both  the  vacuum  chamber 
calibration  curve  and  Andersen's  open-atmosphere  calibration  curve  are 
shown  in  Figure  3-4.  The  curves  have  the  same  general  shape,  but  the 
vacuum  chamber  curve  displays  a  larger  response  at  short  wavelength  than 
does  Andersen's  atmospheric  curve.  The  increased  reponse  can  be  attributed 
to  the  removal  of  Shumann-Runge  oxygen  band  absorption.  The  remainder 
of  the  vacuum  calibration  curve  shows  smaller  response  than  Andersen's, 
and  this  is  believed  to  be  due  to  a  decrease  in  the  irradiance  firom  the 
deuterium  lamp.  There  is  also  a  minor  absorption  peak  near  2740  A  in  the 
vacuum  calibration  curve  which  is  attributed  to  shadowing  by  a  small  particle 
on  the  detector  face. 
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Figure  3-3  Change  in  Detector  Response  due  to  Film 
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Figure  3-4  Comparison  of  Vacuiun  Chamber  and  Open  Atmosphere 
Calibration  Ciurves  for  MUSTANG 
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IV.  ANALYSIS 


A  total  of  8400  spectra  were  recorded  in  two  altitude  groups:  from  150 
km  to  320  km  during  ascent,  and  from  200  km  to  100  km  during  descent.  To 
facilitate  their  evaluation,  the  spectra  were  averaged  into  25  representative 
spectra  each  corresponding  to  a  10  km  segment  of  altitude.  The  averaged 
spectra  are  labelled  with  the  median  altitude  for  the  segment,  e.g.  the  spectra 
from  170  km  to  180  km  are  designated  the  175  km  spectrum.  The  number  of 
individual  spectra  averaged  together  in  an  altitude  segment  varies,  because  of 
the  changing  vertical  speed  of  the  rocket,  from  81  at  105  km  to  330  at  315  km. 
The  resulting  25  spectra  are  herein  modeled  with  synthetic  spectra  from 
which  LBH  column  emission  rate  and  nitric  oxide  density,  can  be  inferred. 

A.  FITTING  THE  OBSERVED  SPECTRA 

During  the  data  collection  period,  a  defective  timing  signal  caused  the 
loss  of  every  seventeenth  data  word  from  the  telemetry  stream.  The  periodic 
loss  of  data  points  caused  a  progressive  shift  in  the  recorded  wavelength 
positions  of  subsequent  points.  The  loss  was  corrected  during  initial  data 
reduction  by  inserting  an  average  value  in  place  of  each  missing  point  and 
shifting  the  positions  of  the  remaining  points.  As  reported  by  Clayton  (1990), 
an  additional  wavelength  discrepancy  caused  the  emissions  from  several 
transitions  to  appear  at  incorrect  wavelengths.  This  discrepemcy  is  thought  to 
be  jm  artifact  of  the  telemetry  circuitry,  and  it  has  been  corrected  by 
performing  a  cubic  spline  fit  between  the  recorded  wavelength  data  and  the 
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theoretical  wavelength  curve  of  MUSTANG  as  determined  by  Aiuiersen 
(1990).  Plots  of  the  measured  peaks  as  a  function  the  spline-fit  wavelength 
data  show  excellent  agreement  with  the  published  values  of  the  peak 
wavelengths. 

Having  established  the  correct  wavelength  distribution  for  the  data, 
synthetic  spectra  were  generated  based  on  temperature  of  the  atmosphere. 

The  temperature  profile  used  for  the  synthetic  spectra  was  determined  by 
Clayton  (1990)  from  his  anal3rsis  of  the  ^(1,0)  emission  of  nitric  oxide.  In  order 
to  accoimt  for  instrument  response  characteristics,  the  synthetic  spectra  were 
convolved  with  the  instrument  slit  function  and  multiplied  by  the 
calibration  parameter  described  in  Chapter  m.  Of  the  512  points  in  each 
spectrum,  the  72  corresponding  to  wavelengths  1848  A  to  2070  A  were  selected 
for  fitting  the  S5mthetic  to  the  observed  spectra.  The  spectrum  in  this 
wavelength  region  is  composed  of  four  emissions  as  shown  in  Figure  4-1. 
Fitting  was  done  using  a  parameter  space  grid  search  adapted  by  Clayton  (1990) 
from  Bevington  (1969).  The  parameters  included  in  the  search  were  LBH 
intensity,  the  sum  of  the  y,  6,  and  e  intensities  from  NO,  and  a  constant 
background  term.  The  background  term  was  included  to  account  for 
scattering  into  the  line  of  sight  from  the  bright  portion  of  the  atmosphere 
below  the  instrument.  The  background  component  is  omitted  from  Figure  4- 
1  for  the  sake  of  clarity.  Goodness  of  fit  was  determined  by  minimizing  the 
value  of  X‘‘  between  the  observed  spectrum  and  the  sum  of  the  S3mthetic 
terms.  The  results  of  the  fitting  are  shown  in  Figures  4-2  to  4-26. 

The  efiect  of  an  exponentially  decreasing  density  is  easily  seen  in  the  fits. 
The  intensity  of  the  emissions  varies  from  nearly  1500  R/A  at  105  km  to  less 
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Comporison  of  Theory  with  Data  at  135  km 
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Figure  4-5  Comparison  of  Synthetic  and  Observed  Spectra  at  135  km 
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Figure  4-6  Comparison  of  Synthetic  and  Observed  Spectra  at  145  km 
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Comparison  of  Theory  with  Data  at  165  km  (ascent) 
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Figure  4-9  Comparison  of  Synthetic  and  Observed  Spectra  at  165  km  (ascent) 
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Comparison  of  Theory  with  Data  at  165  km  (descent) 
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Figure  4-10  Comparison  of  Synthetic  and  Observed  Spectra  at  165  km 
(descent) 
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Figure  4-12  Comparison  of  Synthetic  and  Observed  Spectra  at  175  km 
(descent) 
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Figure  4-13  Comparison  of  Synthetic  and  Observed  Spectra  at  185  km 
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Figure  4-14  Comparison  of  Synthetic  and  Observed  Spectra  at  195  km 
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Figure  4-15  Comparison  of  Synthetic  and  Observed  Spectra  at  205  km 
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Figure  4-17  Comparison  of  Synthetic  and  Observed  Spectra  at  225  km 
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Figure  4-18  Comparison  of  Synthetic  and  Observed  Spectra  at  235  km 
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Figure  4-19  Comparison  of  Synthetic  and  Observed  Spectra  at  245  km 
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Figure  4-23  Comparison  of  Synthetic  and  Observed  Spectra  at  285  km 
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Figure  4-24  Comparison  of  Synthetic  and  Observed  Spectra  at  295  km 
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Figure  4-25  Comparison  of  Synthetic  and  Observed  Spectra  at  305  km 
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Figure  4-26  Comparison  of  Synthetic  and  Observed  Spectra  at  315  km 
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than  12  R/A  at  315  km.  Such  a  large  change  in  the  magnitude  of  the 
measured  quantity  causes  problems  with  signal-to-noise  ratio  as  can  be  seen 
by  comparing  Figures  4-2  and  4-26.  At  105  km,  nearly  every  featvue  of  the 
observed  spectrum  is  accounted  for  in  the  synthetic  spectrum.  Note  that  the 
spectrum  varies  smoothly  and  has  almost  no  spikes.  In  contrast,  the  315  km 
spectrum  shows  rapid  fluctuation  with  wavelength  and  has  many  spikes. 

The  sn\all  magnitude  of  the  observed  intensity  exaggerates  the  effect  of  minor 
variations.  Correspondingly,  there  are  several  features  which  are  not 
mirrored  in  the  s)mthetic  fit.  Nonetheless,  the  general  shape  of  the  315  km 
spectrum  is  well  matched  by  the  sjmthetic  spectrum,  and  prominent 
emissions  such  as  the  LBH  (3,9)  at  1854  A  and  the  y  (2,0)  at  2047  A  are 
accurately  modelled.  The  effect  of  altitude  on  the  relative  intensities  between 
the  LBH  and  the  nitric  oxide  emissions  can  also  be  seen  in  the  figures.  Their 
different  excitation  mechanisms  cause  the  NO  and  N2  emissions  to  dominate 
at  different  altitudes.  The  NO  enussions  are  much  more  intense  than  the 
LBH  enussions  at  low  altitudes.  Emitted  intensities  from  the  two  species  are 
comparable  around  175  km  altitude,  with  the  LBH  contributions  becoming 
dominant  above  185  km  as  photoelectron  flux  increases. 

After  the  ir.Hial  fitting  was  complete,  the  plots  were  evaluated  for  trends 
which  would  indicate  ways  to  improve  the  fits.  In  comparing  the  synthetic 
spectra  to  the  observed,  it  was  noticed  that  the  contributions  from  some 
vibrational  transitions  were  consistently  overestimated  while  other  were 
underestimated.  When  the  intensities  at  a  given  peak  did  not  coincide,  the 
difference  between  the  S3mthetic  and  observed  intensities  as  a  fraction  of  the 
total  intensity  was  used  as  a  correction  to  the  synthetic  intensity.  As  discussed 


67 


in  Chapter  n,  the  influence  of  the  vibrational  terms  on  a  spectrum  is 
determined  by  the  Franck-Condon  factors.  Correcting  the  intensity  value 
generated  by  the  synthetic  spectrum  model  is  equivalent,  then,  to  changing 
the  Franck-Condon  factor  for  the  vibrational  transition  concerned.  Once  the 
correction  factors  were  foimd,  new  S3mthetic  spectra  were  created  and  run 
through  the  fitting  process.  The  specific  vibrational  transitions  used  in  the 
comparisons  are  listed  in  Table  1  along  with  the  scale  factors  they  generated. 
Figure  4-27  gives  a  graphic  representation  of  the  corrections  to  nine  of  the 
LBH  bands  as  a  function  of  vibrational  levels. 

Generation  of  the  correction  factors  was  complicated  by  two 
characteristics  of  the  spectra.  First,  there  were  some  cases,  for  example  the 
e(0,0)  band  at  1877  A,  in  which  the  synthetic  intensity  was  not  consistent  in  its 
deviation  from  the  observed  intensity.  When  this  happened,  the  sjmthetic 
intensity  was  left  unchanged,  or  was  corrected  to  fit  the  largest  possible 
number  of  the  observed  spectra.  Second,  and  more  importantly,  all  four  of 
the  band  systems,  NO  Y/  8/  6/  ai'd  LBH,  exhibit  some  degree  of  overlap  in  the 
wavelength  range  of  interest.  This  overlap  means  that  the  corrections  for  the 
bands  are  coupled.  In  an  effort  to  reduce  the  effect  of  the  coupling,  peaks  in 
the  synthetic  NO  spectrum  were  compared  with  the  measured  spectrum  at 
105  km,  where  the  nitric  oxide  emissions  are  much  greater  than  those  from 
N2.  Similarly,  the  LBH  comparisons  were  made  with  the  265  km  spectrum. 
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TABLE  1 


CORRECTIONS  TO  FRANCK-CONDON  FACTORS 

Vibrational  Bandhead 

Band  Correction  Wavelength  (A) 

N2  LBH 


1,8 

1.3 

1895 

2,9 

1.0 

1914 

2,10 

1.0 

1992 

3,9 

1.2 

1854 

3,10 

0.8 

1929 

3,11 

0.8 

2008 

4,10 

1.0 

1873 

4,11 

0.9 

1948 

4,12 

0.9 

2026 

5,11 

1.4 

1889 

5,12 

0.75 

1964 

5,13 

1.2 

2045 

6,13 

1.0 

1983 

2,0 

1.0 

2047 

3,0 

1.0 

1956 

NO  6 

No  adjustments  made  to  individual  Franck-Condon  factors. 

NOe 


0,0 

0.9 

1877 

0,1 

1.15 

1945 

1,1 

0.6 

1863 
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Figure  4-27  Corrections  to  Franck-Condon  Factors  for  Selected  LBH 
Transitions 
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for  LBH  (v*,v")  transition 


Figure  4>27  indicates  that  the  Franck-Condon  factors  measured  in  this 
experiment  agree  fairly  closely  with  theory  for  the  (4,v")  transitions.  It  also 
shows  that  the  (3,9),  (5,11),  and  (5,13)  transitions  are  underestimated  by  theory, 
while  the  (3,10),  (3,11),  and  (5,12)  transitions  are  overestimated.  These 
disagreements  between  measured  and  theoretical  values  may  be  partly  due  to 
the  influence  of  strong  emissions  from  the  nitric  oxide  e(l,l)  and  e(0,0)  bands 
which  contaminate  the  portion  of  the  spectrum  from  which  the  LBH 
corrections  were  estimated. 

B.  COLUMN  EMISSION  RATES 

When  the  fits  were  satisfactory,  the  scale  factors  between  synthetic  and 
observed  spectra  represented  the  colunm  emission  rates  for  nitric  oxide  and 
nitrogen.  This  rate  can  be  converted  directly  to  a  column  density  for  NO.  The 
density  profile  for  NO  is  shown  in  Figure  4-28.  The  error  bars  represent 
estimates  of  the  maximum  and  minimum  scale  factors  which  will  produce  a 
plausible  fit  between  the  model  and  the  observed  data.  The  profile  shows  the 
expected  exponential  decrease  of  density  with  increasing  altitude. 

The  results  from  the  model  calculations  of  LBH  emissions  are  shown  in 
Figure  4-29.  These  results  represent  the  column  emission  rate  profile  for  all 
of  the  LBH  bands  of  N2.  These  bands  extend  over  a  wavelength  range  of  1400 
A  from  1205  A  to  2602  A.  In  order  to  compare  the  emission  rates  from  the 
MUSTANG  data  to  the  model,  it  was  necessary  to  calculate  the  amount 
contributed  to  the  entire  LBH  band  system  emission  by  the  thirteen  bands  in 
the  1850  A  to  2070  A  region.  This  calculation  was  performed  by  summing  the 
Fremck-Condon  factors  for  the  thirteen  bands  listed  in  Table  1  and  dividing  by 
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Figxire  4-28  Nitric  Oxide  Column  Density 
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the  sum  of  all  Franck-Condon  factors  for  the  LBH  S3rstem.  These  bands 
represent  approximately  21.1  percent  of  the  total  LBH  emissions.  The  error 
estimates  for  the  LBH  emission  rates  again  represent  the  maximum  and 
minimum  scale  factors  for  achieving  a  believable  fit  with  the  observed 
spectra.  The  estimated  error  for  the  LBH  terms  is  relatively  large  at  low 
altitude  because  the  LBH  contribution  to  the  recorded  spectrum  is 
overwhelmed  by  the  nitric  oxide  contributions.  At  high  altitude,  the  small 
signal-to-noise  ratio  leads  to  relatively  large  error  estimates  again. 

As  noted  in  Chapter  n,  the  column  emission  rate  for  LBH  cannot  be 
converted  directly  to  an  N2  density  because  the  photoelectron  flux  depends  on 
that  density.  Consequently,  Equation  2.15  caimot  be  solved  directly  for  the 
density  of  N2.  Instead,  the  forward  calculation  must  be  iterated  by  varying  the 
nitrogen  density  and  the  photoelectron  flux  until  the  calculated  emission  rate 
agrees  with  the  observed  rate.  The  photoelectron  flux  was  estimated  using 
two  computer  programs  for  atmospheric  modelling.  The  neutral  atmosphere 
was  modelled  by  MSIS-83  (Hedin,  1983)  using  environmental  conditions 
pertinent  to  the  MUSTANG  launch  as  inputs.  The  ouput  from  MSIS  is  a 
model  of  the  neutral  atmosphere  with  values  for  N2  density,  O2  density,  O 
density,  and  atmospheric  temperature  as  a  function  of  altitude.  The  model 
atmosphere  was  then  fed  into  the  NRL  PEGFAC  model  which  calculated  the 
photoelectron  flux  using  a  Monte  Carlo  techniques.  The  diamond  symbols  in 
Figure  4-29  indicate  the  column  emission  rate  measured  by  MUSTANG.  The 
MUSTANG  curve  clearly  does  not  coincide  with  the  model  curve.  There  are 
a  number  of  explanations  for  the  disparity  between  the  profiles,  one  of  which 
may  be  a  difference  between  the  actual  temperahure  profile  and  that  used  in 
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the  models.  The  temperature  at  high  altitude  in  the  thermosphere  is  called 
T-  infinity  or  Tinf.  T^f  calculated  by  MSIS  was  1312  K.  Clayton  (1990),  in  his 
evaluation  of  the  MUSTANG  data,  found  a  value  for  Tinf  1800  K.  The 
value  of  Tinf  ^  large  impact  on  calculated  column  emission  rate,  especially 
at  high  altitude.  The  effect  is  shown  in  Figure  4-30  which  presents  emission 
rates  calculated  at  1312  K,  15(X)  K,  and  1700  K.  The  difference  in  the  two 
values  may  be  partly  explained  by  the  differing  methods  used  to  calculate  Tinf. 
Clayton  produced  his  temperature  profile  by  comparing  the  bandwidth  of  the 
NO  y(1/0)  emission  from  a  synthetic  spectrum  with  that  of  the  observed  y(1,0) 
emission.  The  MSIS  model  uses  the  FI  0.7  solar  activity  indication  to  create  its 
temperature  profile.  The  F10.7  value  is  a  measure  of  solar  output  in  the  radio 
frequency  portion  of  the  electromagnetic  spectrum  which  has  become  a 
benchmark  for  measuring  solar  activity.  Unfortunately,  changes  in  F10.7  are 
not  linearly  proportional  to  changes  in  UV  activity,  so  the  F10.7  index  may 
not  indicate  the  same  temperature  profile  as  that  indicated  by  analysis  of  the 
bandwidth  of  ultraviolet  emissions.  Because  his  temperature  profile  was 
based  on  analysis  of  the  MUSTANG  data,  Clayton's  value  of  Tf^f  was  used  in 
this  analysis. 

Another  reason  for  the  disparity  is  the  influence  of  the  profiles  for 
nitrogen  density,  molecular  oxygen  density,  and  atomic  oxygen  density, 
written  [N2I,  [O2I/  and  [O]  respectively.  The  MSIS  density  profiles  zu-e  based  on 
multi-year  observations  of  densities  under  conditions  similar  to  those  given 
as  inputs  to  the  model.  While  this  approach  would  be  expected  to  produce 
reliable  average  results,  individual  cases  may  vary  significantly  fi'om  the  long 
term  average.  The  [N2]  and  [O]  profiles  used  in  the  PEGFAC  models  have  a 
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Figure  4>30  Effect  on  PEGFAODerived  Emission  Profile  of  Variation  in  Tint 
Parameter  of  MSIS 
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great  influence  on  the  calculated  column  emission  rate  as  can  be  seen  in 
Figures  4-31  and  4-32.  The  solid  curve  in  each  figure  represents  the  emission 
rate  derived  from  PEGFAC  using  the  density  calculated  by  MSIS,  while  the 
other  curves  show  the  effect  of  a  30%  variation  from  the  MSIS  value.  Figure 
4-31  indicates  that  a  change  of  30%  in  the  density  of  nitrogen  at  120  km  causes 
the  emission  rate  to  vary  by  approximately  30%  at  high  altitudes.  The  peak 
value  of  intens‘*y  varies  by  only  about  5%,  but  the  height  of  the  peak 
increases  with  increasing  nitrogen  density.  The  oxygen  density  profile  in 
Figure  4-32  shows  quite  different  properties.  A  30%  change  in  [O]  has 
negligible  effect  on  emission  rate  above  250  km,  but  causes  the  peak  emission 
rate  to  vary  by  nearly  20%.  The  peak  value  decreases  with  increasing  oxygen 
density,  and  the  altitude  of  the  peak  remains  constant.  Changing  the  density 
of  O2  has  little  effect  on  the  emission  rate  above  110  km  (Figure  4-33). 

Two  different  emission  rate  profiles  from  PEGFAC  are  compared  with 
the  MUSTANG  profile  in  Figures  4-34  and  4-35.  The  MSIS  parameters  used 
to  generate  Figure  4-34  are:  Tinf  =  1500  K,  densities  at  120  km;  [N2]  =  4.10  x  lO'^ 
cm'^,  [O]  =  6.36  x  10^®  cm*^,  [O2]  =  4.66  x  10^®  cm"^.  In  Figure  4-35  the  values 
are:  T^f  =  1700  K,  [N2]  =  3.15  x  10^^  cm"^,  [O]  =  4.55  x  10^®  cm"^,  and  [O2]  =  4.66  x 
10^®  cm‘^.  While  both  figures  indicate  fairly  good  agreement  in  overall  shape 
between  the  curves,  in  neither  case  zue  all  features  of  the  MUSTANG  profile 
accounted  for  by  the  model.  The  peak  of  the  intensity  profile  is  not  well  fitted 
in  both  cases.  Above  195  km,  the  model  in  Figure  4-34  shows  good  agreement 
with  the  MUSTANG  data.  Unfortunately,  this  agreement  indicates  a  value 
for  Tinf  of  1500K,  whereas  Clayton's  analysis  cf  band  widths  indicated  Tinf  is 
1800  K.  Figure  4-35,  at  1700  K,  shows  that  the  fit  at  high  altitude  is  worse  as 
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Tinf  approaches  Clayton's  value  of  1800  K.  In  addition,  the  values  of  [O]  in  the 
two  figures  are  not  based  on  an  independent  measure  of  the  atomic  oxygen 
density.  When  it  was  clear  that  the  MSIS  values  for  oxygen  and  nitrogen 
densities  would  not  give  a  close  fit  to  the  data,  density  vsilues  were  chosen 
arbitrarily  in  an  attempt  to  fit  the  model  profile  to  the  MUSTANG  profile. 

The  fact  that  Figure  4-35  shows  reasonable  agreement  with  the  observed 
profile  despite  the  large  differences  in  parameter  values  from  those  in  Figure 
4-34  indicates  that  the  method  used  in  this  thesis  will  not,  by  itself,  allow 
calculation  of  a  unique  value  for  nitrogen  density. 
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Figiire  4-31  Effect  on  PEGFAC-Derived  Emission  Rate  of  Variation  in  N2 
Density  Parameter  of  MSIS 
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Figure  4-32  Effect  on  PEGFAC-Derived  Emission  Rate  of  Variation  in  O 
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Figure  4-34  Comparison  of  Coliunn  Emission  Rate  Profiles.  The  MSIS  curve 
is  for  the  following  conditions:  E>ensities  at  120  km;  IN2]  =  4.10  x  10^^  cm'^, 
[OJ  =  6.36  X  10^°  cm-3  [02!  =  4.66  x  10*0  cm'^;  =  1500  K. 
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Figure  4-35  Comparison  of  Column  Emission  Rate  Profiles.  The  MSIS  curve 
is  for  the  following  conditions:  Densities  at  120  km;  [N2]  =  3.15  x  10^^  cm'3, 
[01  =  4.55  X 10^°  cm-3  [O2]  =  4.66  x  10^0  cm-3;  Tmf  =  1700  K. 
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V.  CONCLUSION 


A.  SUMMARY  OF  RESULTS 

This  thesis  has  demonstrated  that  synthetic  spectra  based  on  the 
calculation  of  excited  state  populations  of  molecular  nitrogen  and  nitric  oxide 
can  be  made  to  model  accurately  the  observed  airglow  in  the  1850  A  to  2070  A 
wavelength  region.  The  scale  factors  between  the  models  and  the  measured 
spectra  produce  reasonable  values  for  line-of-sight  column  densities  of  nitric 
oxide.  Column  densities  for  molecular  nitrogen  caimot  be  determined 
uniquely  using  the  method  outlined  in  this  thesis.  A  vacuum  calibration 
using  a  chamber  manufacttued  from  PVC  pipe  has  been  presented.  The 
calibration  was  performed  in  a  nitrogen  atmosphere  to  prevent  absorption  of 
the  ultraviolet  radiation  by  oxygen.  Cham^  »r  geometry  and  hardware 
considerations  dictated  that  the  vacuum  calibration  curve  be  used  to  extend 
an  existing  absolute  calibration  curve. 

B.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

The  utility  of  measiuring  LBH  emissions  as  a  means  to  find  nitrogen 
density  depends  critic2dly  on  the  ability  to  model  the  photoelectron  flux 
which  drives  those  emissions.  An  independent  measure  of  the  density  of 
atomic  oxygen  is  essential  to  reducing  the  number  of  free  parameters  used  in 
generating  the  theoretical  coliunn  emission  rate  for  the  LBH  bands.  An 
acoirate  temperature  profile  of  the  upper  atmosphere  is  also  needed.  Further 
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insight  into  trends  in  the  corrections  to  Franck-Condon  factors  may  be  gained 
by  fitting  S3mthetic  spectra  over  a  wider  wavelength  range  the  has  been  done 
to  date. 

Vacuum  calibration  of  the  MUSTANG  spectrograph  can  be  improved  by 
repeating  the  procedure  in  a  stainless  steel  vacuum  chamber.  A  stainless  steel 
chamber  would  eliminate  the  problem  of  outgassing  experienced  with  PVC, 
and  would  allow  the  recording  of  spectra  under  high  vacuum. 
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